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We derive the generalized Markovian description for the non-equilibrium Brownian motion of a 
heated particle in a simple solvent with a temperature-dependent viscosity. Our analytical results 
for the generalized fluctuation-dissipation and Stokes-Einstein relations compare favorably with 
measurements of laser-heated gold nano-particles and provide a practical rational basis for emerging 
photothermal technologies. 
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Brownian motion is the erratic motion of suspended 
particles that are large enough to admit some hydrody- 
namic coarse-graining, yet small enough to exhibit sub- 
stantial thermal fluctuations. Such mesoscopic dynamics 
is ubiquitous in the micro- and nano- world, and in partic- 
ular in soft and biological matter [TJ |2] . Since their first 
formulation more than a century ago, the laws of Brow- 
nian motion have therefore found so many applications 
and generalizations in all quantitative sciences that one 
may justly speak of a "slow revolution" 3J. In Langevin's 
popular formulation they take the simple form of New- 
ton's equation of motion for a particle of mass m and 
radius R subject to a drag force —^oP/ m and a randomly 
fluctuating thermal force £(t); 



p + CoP/m = £ (t > 0) 



(1) 



As a cumulative representation of a large number of 
chaotic molecular collisions £ is naturally idealized as a 
Gaussian random variable. Its variance is tied to the 
Stokes friction coefficient 



Co = Qtt7] R 



(2) 



in a solvent of viscosity 770 such as to guarantee consis- 
tency of the averages (. . .) over force histories with 
Gibbs' canonical ensemble, namely 

(£(*)) = , (Ci(t)C-(O)) - 2* fl T CoMW ■ ( 3 ) 

This prescription implements the fluctuation-dissipation 
theorem for the system comprising the Brownian parti- 
cle and its solvent at temperature T . Strictly speaking, 
in view of how it deals with long-ranged and long-lived 
correlations arising from conservation laws governing the 
solvent hydrodynamics, this practical and commonplace 
Markovian description applies only asymptotically for 
late times [H[S]. Corresponding corrections to Eqs. (2][3) 
are accessible to modern single-particle techniques and 
become most relevant in nano-structured environments 

®G|. 

Thanks to its prominent role in the "middle world" [2] 
between macro- and micro-cosmos, and its experimen- 
tal and theoretical controllability, Brownian motion has 



become a "drosophila" for formulating and testing new 
(and sometimes controversial) developments in equilib- 
rium and non-equilibrium statistical mechanics [51414). 
In this Letter, we introduce a non-equilibrium general- 
ization that has so far received little attention, namely 
the Brownian motion of a particle maintained at an el- 
evated temperature T p > T . From its hypothetical 
sibling ("cool Brownian motion", T p < Tq) such "hot 
Brownian motion" (HBM) is distinguished by having ob- 
vious realizations of major technological relevance such 
as nano-particles suspended in water and diffusing in a 
laser focus. Due to a time-scale separation between heat 
conduction and Brownian motion these particles carry 
with them a radially symmetric hot halo easily detected 
with a second laser. This provides the basis for promis- 
ing photothermal particle tracking |15j and correlation 
spectroscopy ("PhoCS") [ToTJIB] techniques with a high 
potential of complementing corresponding fluorescence 
techniques 19 in numerous applications. However, a 
photothermal measurement necessarily disturbs the dy- 
namics it aims to detect more severely than typical flu- 
orescence measurements, so that the development of an 
accurate theoretical description of the Brownian motion 
of heated particles is a crucial prerequisite for making the 
method competitive. This is not an entirely straightfor- 
ward task (as some might suggest |30p and requires an 
extension of the familiar theory, as explained in the fol- 
lowing. We arrive at simple analytical generalizations of 
Eqs. (2][3), which should be sufficiently accurate for most 
practical applications. 

For clarity, we restrict the following discussion to an 
idealized situation: a hot spherical Brownian particle of 
radius R at the center of a co-moving coordinate sys- 
tem in a solvent with a temperature-dependent viscosity 
r](T) that attains the value rjo at the ambient tempera- 
ture Tq imposed at infinity. Favorable conditions are as- 
sumed, such that potential complications resulting from 
long-time tails [7], convection [2U], thermophoresis [3T], 
etc. can be neglected. To avoid confusion in comparisons 
with experimental data, we do however distinguish the 
solvent temperature T s at the hydrodynamic boundary 
corresponding to the particle surface from the particle 
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temperature T p itself, as these may differ substantially 
[22] . It is the temperature difference AT = T s — To 
that determines the heat flux responsible for the non- 
equilibrium character of the problem. On relevant time 
scales, the resulting temperature field around the particle 
follows from the stationary heat equation, i. e. 

T(r) = T + BAT/r . (4) 

The task of finding appropriate generalizations of Eqs. ^ 
[3]) under these conditions is split into two steps corre- 
sponding to the two force terms in Eq. ([I]), the damping 
and the driving force, or friction and thermal noise, re- 
spectively. 

The first goal is mainly technical, namely to generalize 
Eq. |2| by solving 

Vp = V • j)(r)[Vu + (Vu) T ; 







(5) 



for the stationary fluid velocity field u(r) under the usual 
no-slip boundary condition. The new feature compared 
to Stokes' classical derivation is the radially varying vis- 
cosity rj(r) resulting from Eq. Q. A numerically precise 
solution of Eq. ^ can be obtained with a differential 
shell method [53] along the lines of similar work for inho- 
mogeneous elastic media [24] . However, for our present 
purposes, as well as for practical applications, we wish to 
find a generally applicable analytically tractable approx- 
imation. We therefore resort to a toy model that evades 
the technical difficulties related to the vector character 
of the fluid velocity but retains the long-ranged nature 
of the hydrodynamic flow field. We replace u(r) by a 
fictitious diffusing scalar u(r) without direct physical sig- 
nificance, for which Eq. ([5| is readily solved analytically. 
More explicitly, Eq. ^ reduces to V • rj(r)X7u(r) = in 
the scalar model. A separation ansatz u(r) = u r (r)u^(i}) 
leads to the radial equation 



[d r + 2/r + (d r ln?7)]<9 r u r = 



(6) 



solved by d r u r cx (ryr 2 ) -1 for physically reasonable func- 
tions rj(r). The quantities u r and r\d r u r are now in- 
terpreted as the analogue of the velocity of the particle 
and the hydrodynamic drag force per area, respectively. 
The generalized effective friction coefficient Chbm of hot 
Brownian motion is then estimated up to a numerical fac- 
tor as their ratio, disregarding the contribution from the 
angular part. A comparison with Eq. ^ in the isother- 
mal limiting case of constant viscosity rj(r) = 770 helps to 
calibrate the model and fix the undetermined numerical 
factor, which is then taken over to situations with radi- 
ally varying rj{r). The accuracy of this procedure can 
be assessed and further improved by a comparison with 
analytical and numerical results from the mentioned dif- 
ferential shell method [23]. Some technical details are 
provided in [25] and the result is summarized in Fig. [I] 

An analytically tractable expression for the effective 
friction coefficient Chbm as a function of temperature fi- 
nally results from a combination of the calibrated model 



▲ >/„[l-(l-l/K)exp(l-r)] 
■ %[l~(l~l/*)/r] 

Vogel-Fulcher, Eq. (7) 
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FIG. 1: Comparison of analytical predictions of the scalar toy 
model to numerical results from the differential shell method 
(symbols) for exemplary long-ranged radial viscosity profiles 
■q(r) with k = r]{r — s> <x)/rj(r = R) in a parameter regime of 
potential practical interest. While the simplest version of the 
model (dashed lines), which employs a constant calibration 
factor, exhibits noticeable systematic errors, the more elabo- 
rate version (solid line), corresponding to Eq. ( 12 1, should be 
sufficiently accurate for practical applications. 



with Eq. Q and a phenomenological expression for the 
temperature dependence of the solvent viscosity such as 



j?(T) = 7?oo exp[A/(T - T VF ) 



(e.g. for water; but power-laws could be processed just as 
well) . The effective friction can be reinterpreted in terms 
of an effective solvent viscosity ?7hbm = CHBM/67ri? that 
replaces r/o in Eq. ([2]) under non-isothermal conditions. 
For reduced temperature increments 9 = AT/(T — 
Tvf) < 1 the result is well approximated by its trun- 
cated Taylor series [2"5] 
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(8) 



This provides the wanted generalization of Eqs. (T][2). 

To turn Eqs. ([l]|3| into a fully predictive Markov model 
of hot Brownian motion, the remaining task is to com- 
pute, in the same spirit, an appropriate effective temper- 
ature Thbm to replace To in Eq. In other words, 
we aim at establishing a generalized non- equilibrium 
fluctuation-dissipation relation for Brownian motion in a 
co-moving radial temperature gradient. In analogy to the 
better understood situation in globally isothermal non- 
equilibrium steady states [26], we expect to retrieve the 
fluctuation-dissipation relation only after excluding the 
"housekeeping heat" from the entropy balance; i.e. the 
heat constantly flowing from the particle to infinity to 
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maintain the temperature gradient. All we have to con- 
sider is the minuscule excess dissipation associated with 
the damped motion of the Brownian particle. In this 
respect, it is crucial to appreciate the long-range corre- 
lated character of the hydrodynamic flow, which affects 
both dissipation and thermal fluctuations. It also helps 
in setting up a systematic coarse-grained calculation by 
extending the standard framework of fluctuating hydro- 
dynamics [27] to moderate temperature gradients |23j . 

In simple terms, the process of Brownian motion can be 
rephrased as a constant transformation of some thermal 
energy from the solvent into an equal amount of kinetic 
energy for the Brownian particle and vice versa. In a 
stationary situation the mutual energy transfer must be 
balanced to obey the first law. More precisely, the spatial 
integral over the local excess dissipation q(r) — i. e. the 
heat created (per unit of time) by the solvent flow at 
position r in response to the movement of the Brownian 
particle — must on average match the rate of kinetic 
energy transfer W p to the particle, 



(W p ) = / dr (q(r)) 



(9) 



Moreover, to respect the second law, the motion must not 
cause a net average entropy change, which was the origin 
of major reservations against the modern interpretation 
of Brownian motion till the early 20th century. However, 
this only means that one has to make sure that the in- 
tegral over the local entropy flux to the solvent — i. e. 
the local dissipation rate q(r) divided by the local solvent 
temperature T(r) — equals on average the entropy flux 
S p = Wp/TuBM conferred to the Brownian particle: 



dr 



<9(r)> = 1 
T(r) Thbm 



dr (q(r)) 



(10) 



This then defines the wanted effective Brownian temper- 
ature Thbm, if the dissipation q(r) is expressed in terms 
of the local viscosity ?y(r) and Vu(r). Within our scalar 
model q(r) — r)(r)[d r u r (r)] 2 /2, hence 



T 



HBM 



dr V (r)((d r u r ) 2 ) 



dv^\{(d r u r f) . (11) 
T(r) 



For the special case of a temperature-independent con- 
stant viscosity rj this reduces to the simple explicit ex- 
pression 



Thbm = AT/ ln(l + AT/r ) 



(12) 



The analytical expression generalizing this to the main 
case of interest, a viscosity r](r) that varies radially ac- 
cording to Eqs. @ & 0, is given in Ref. |25j . For small 
temperature increments AT <C Tq a practical approxi- 
mation is 

Thbm « T +AT/2-[l - Hvo/Voo)} AT 2 /(24T ) . (13) 
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FIG. 2: The effective diffusion coefficient D H bm(AT) of hot 
gold nano-particles traversing a laser focus in water: experi- 
mental data (open/closed symbols for R — 40/60 nm) versus 
analytical predictions from the scalar model (solid lines); for 
solvent and focus parameters and the error bars see Ref. [25]. 
Inset: Dhbm(AT) according to numerical predictions from 
the differential shell method (squares), analytical solutions of 
the scalar model (corresponding to the lowest pair of curves 
in Fig. and the naive suggestion to identify the HBM pa- 
rameters with the conditions at the particle surface (dotted); 
the agreement between the symbols and the solid line demon- 
strates the equivalence of Eqs. ( 10 1 & (111. 



A hot Brownian particle described by Eqs. (|2j3| with 
rjQ and Tq replaced by the corresponding effective quan- 



tities ijhbm and Thbm from Eqs. (12) & (15) performs a 



random diffusive motion characterized by an effective dif- 
fusion coefficient Z?hbm obeying the generalized Stokes- 
Einstein relation 



D 
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fesTlBM 
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(14) 



This prediction is tested against the numerical differen- 
tial shell method in the inset of Fig. [2j The good agree- 
ment demonstrates the equivalence of Eqs. ( 10 1 and ( 11 ) 



In order to test Eq. ( 14 ) also experimentally, we used a 



photothermal microscopy setup with gold nano-particles 
in water, as described in Refs. [THIUS]. Particles passing 
through the common focal volume of a heating and a de- 
tection laser beam leave a trace of photothermal bursts 
in the detector, which encodes information about the dif- 
fusivity. The spatially inhomogeneous heating power in 
the laser focus implies, via Eq. fl4| ), that the diffusion 
in the focus is inhomogeneous. (A wider focus of the 
heating laser would avoid this complication but is gener- 
ally undesirable as one wants to minimize sample irradia- 
tion.) We therefore pursue a first-passage time approach 
to determine the apparent effective diffusion coefficient 
T*hbm °f inhomogeneous hot Brownian motion from the 
burst durations, which we identify with the transit times 
of the particles passing through the focus volume [31] . 



4 



The time periods r during which the photo-thermal sig- 
nal supersedes a fixed percentage of the maximum signal 
at a given laser power are recorded for a large number of 
photothermal bursts. The diffusion coefficient is then ex- 
tracted from the exponential decay of the obtained tran- 
sit time distribution P(r) at large r [2"51 l2"gl |2"5 ] . 



In P(t -> oo) oc — -Dhbm 7 " • 



(15) 



Figure [2] shows the result of such measurements for vari- 
ous laser powers. The surface temperatures T s = T +AT 
have been calculated from known quantities, namely the 
incident laser intensity, the optical absorption coefficient 
of the particles, and the heat conductivity of the solvent 
[18] . Due to our limited knowledge of the focus geometry, 



the factor of proportionality in Eq. ( 15 1 could not be de- 



termined precisely, though. We therefore took the liberty 
to multiply each data set by an overall factor to optimize 
the fit [25]. Yet, the good agreement of the functional de- 
pendence with the prediction provides strong support for 
our analytical results, over a considerable temperature 
range. At the same time, it establishes hot Brownian 
motion as a robust and manageable tracer technique. 

In summary, by introducing appropriate effective 
friction (viscosity) and temperature parameters Chbm 
(?7hbm) an d Thbmi f° r w hi cri we provided explicit an- 



alytical expressions in Eqs. (12 1 and (15), the convenient 



Markovian description of Brownian motion in terms of 
Eqs. (2|3| could be extended to non-equilibrium condi- 
tions, where the temperature of the Brownian particle 
differs from that of the solvent. While Eqs. ([2J3]) are re- 
covered in the isothermal limit, the general predictions 
differ significantly from what might have been guessed 
from simple rules of thumb and provide an instructive 
illustration of the general dictum that hydrodynamic 
boundary conditions should not be confused with the mi- 
croscopic conditions at the boundary [8] . We sidestepped 
some technical difficulties of the corresponding problem 
in fluctuating hydrodynamics by introducing an analyti- 
cal toy model that we calibrated with help of more elab- 
orate analytical and numerical calculations. Our analyti- 
cal prediction for the effective diffusion coefficient, based 



on the generalized Stokes-Einstein relation in Eq. (14), 



compares favorably with our measurements of gold nano- 
particles depicted in Fig. [2] and thus provides a conve- 
nient basis for photothermal tracer techniques [T5j [18] 
with a high potential of complementing corresponding 
fluorescence-based methods applied in many fields from 
nano-technology to biology. 

We acknowledge helpful discussions with A. Wiirger, 
and financial support from the Deutsche Forschungsge- 
meinschaft (DFG) via FOR 877. 
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Supplementary Material 

The supplementary material is organized as follows. In Section 1, we revisit Stokes' problem of the viscous drag 
on a sphere for an inhomogeneous viscosity rj{r) approximated by (i) a step function and (ii) a staircase, from 
which we obtain numerically exact solutions for ?7hbm i n the continuum limit. In Section 2, we solve the scalar toy 
model. Comparison with (i) suggests an improved calibration. Section 3 provides the complete expression for the 
diffusion coefficient Z?hbm and Section 4 the generalization -Dhbm for inhomogeneous hot Brownian motion. Section 5 
summarizes some phenomenological parameter values. For complete derivations and a more comprehensive discussion 
see [51]. 



1. SHELL AND DIFFERENTIAL SHELL METHOD FOR STOKES' PROBLEM 

Stokes' classical problem of finding the friction coefficient C of sphere of radius R in a homogeneous Newtonian fluid 
of known viscosity 770 has a well-known solution: £ = 6^770 R. As stated in the main text, we wish to generalize this 
result to radially varying viscosities rj(r). We consider (i) a step profile and (ii) a staircase profile. In both cases the 
solution for the velocity u and pressure p take on the form 

( a i 2 a 3 A 

u r = I a H h a 2 r H — ? cos 

V r r J / 

ug =- (a + ^ + 2a 2 r 2 - ^)sm9 (SI) 
V 2r 2r A J 

P = Po + + 1 0a 2 rj T) cos 9 

in each of the spatially homogeneous regions. 

(i) The step profile r)(r) = r)$ [! + («— l)0(r — b)] jumping from rj s to rjo = Krj s at r = b. The coefficients ao . . . 03 
in the two domains r < b and r > b follow from the continuity conditions for the velocity and the stress at the 
jump discontinuity. 

(ii) Similarly, the solution for the staircase profile is characterized by a set (ao, a±, a 2 , a 3 )j of coefficients and corre- 
sponding continuity conditions. Along the lines of similar work for inhomogeneous elastic media |S2j . we take the 
limit of an infinite staircase of infinitesimally thin shells, and the coefficients become radial functions flj(r) that 
can generally only be evaluated numerically. This differential shell method yields the (inverse) effective friction 
coefficient 

C 1 = a (R) - a^RjR- 1 + a 2 (R)R 2 + a 3 (R)R- 3 \ F=1 , (S2) 

evaluated at a reduced force of F — 1. (The friction coefficient is defined by the force, viz. the stress integrated 
over the particle surface, divided by the velocity relative to the fluid at infinity.) Results for some exemplary 
viscosity profiles are shown as symbols in Fig. I of the main text. 



2. SCALAR TOY MODEL 



To find a generally applicable analytically tractable approximation for the friction coefficient Chbm we resort to a 
toy model that evades the technical difficulties related to the vector character of the fluid velocity but retains the 
long-ranged nature of the hydrodynamic flow field. We replace u(r) by a fictitious diffusing scalar u{r) without direct 
physical significance, for which Eq. (5) reduces to V • 7y(r)Vw(r) = 0. Hence, we seek a radially symmetric solution 
u r of Eq. (6), [d r + 2/r + (d r \mf)]d r u r — 0. Integrating twice, 

f°° dr' 

u r (r) = K K = const. , (S3) 

J r r l r\(r') 

which simplifies to K/rj^r for homogeneous viscosity r)(r) = rio, and can be expressed in the following closed form for 
T(r) = Tq + RAT/r, AT > and the Vogel-Fulcher temperature-dependence of the viscosity specified in Eq. (7) of 
the main text: 

a 

(S4) 

l + fi/r- 



u r (r) 



Ka 



Ei(-x) 
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The abbreviations a = B/Tq and 



RAT/Tq have been used. The effective friction coefficient is 
4TrR 2 T 1 d r u r (R) _ 4nK 



HBM 



U r (R) 



U r (R) 



(S5) 



For homogeneous viscosity rj{r) = ?7o it degenerates to AirrjoR, indicating a mismatch by a factor of 3/2 compared 
to exact result. The simplest calibration of the scalar model consists in correcting this constant factor such as to 
match the predictions in the homogeneous case. In the inhomogeneous case, AT > 0, we express Eq. (^5j) by the 
corresponding effective viscosity ?7hbm- Including the mentioned factor of 3/2 and making the temperature dependence 
(see Section 5) explicit, we have (with the abbreviation T* = T — Tvf) 
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(S6) 



This result features as the lowest dashed line in Fig. 1 of the main text; analogous calculations for different viscosity pro- 
files rj(r) provide the other dashed lines in the figure. Viscosity profiles of the form rj(r) = t]q[1 — (1 — l/n)(R/r) n ], n E IN 
result in 



no 



?7hbm V n 



1 



1 K — 1 

n ' k 



(S7) 



where 2-fi(a, b; c; z) denotes the hypergeometric function. For the two cases n = 1,2, shown in Fig. 1 of the main 
text, 77o/?7hbm can be expressed as 



>h> 



k In k 



?7hbm k 



1 



(n=l) 



(S8) 



and 



»7o 



k arctanhi 



^HBM 



y/ K ( K - x ) 



(n = 2) 



(S9) 



For a more sophisticated calibration of the scalar model, we consider again the step profile (i) . The scalar model 
with the simple calibration predicts the friction coefficient 



c 



b/R 



1 + (&/#-!)« 



67r?7oi? 



(S10) 



By k = r/o/ris, we denote the ratio of the ambient viscosity and the solvent viscosity at the surface of the Brownian 
particle, as before. Note that the trivial limits K — > 1 and b — > R, and the limit k — > of a frozen surface layer, are 
correctly obtained, whereas the joint limit k — > oo and b — > R, corresponding to a particle coated with an infinitesimal 
superfluid layer, is ambiguous. To recover the correct slip boundary condition (£ = Anr/R) in this case, one has to 
take this limit along the curve defined by b/R = 1 + 1/2k. If we impose this constraint on Eq. (£ 10), the calibration 
factor to match it with the exact solution, which we obtain as described in the previous section, viz. Eq. (S^]), is found 
as 



3(80/- 



80k 3 + 60k 2 + 20k + 3) 



200k 4 + 250k 3 + 205k 2 + 65k + 9 



(Sll) 



With this more elaborate calibration, the analytical predictions of the scalar model (solid lines in Figures 1 and 2 of 
the main text) practically coincide quite universally over a broad range of k with the numerical predictions obtained 
from the differential shell method. 

For moderate temperature increments AT «0... 150 K, which are probably of greatest interest in practical appli- 
cations, the result can be further simplified by expanding Eqs. (S^jE 11 ) in a series in 9 = AT/(Tq — Tvf), 
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(S12) 



Truncation after the second order yields Eq. (8) of the main text. In Fig. SI it is compared to the full expression 



and to the result obtained with the simple calibration by the constant factor 3/2. 
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FIG. SI: Temperature dependence of the effective viscosity according to the scalar model, Eq. (q5J): simplest calibration by a 
constant factor to match the isothermal limit (dashed); improved calibration with Eqs. (i^i 11 1 (Sash-dotted); the truncated 
Taylor series of Eq.(tfl2| (solid). 



3. EFFECTIVE DIFFUSION COEFFICIENT 



Inserting the solution u{r) of the scalar model from Eq. (£j4j) into Eq. (11) of the main text yields the effective 
temperature 
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with the following abbreviations 



(S13) 
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(S14) 



Again, a simpler approximate expression is expected to suffice for most practical purposes, in particular with regard to 
the various minor contributions that were neglected altogether in our approach. For moderate temperature increments 
AT <C Tq a Taylor expansion yields 



Thbm 
T n 



= 1 



AT 

Wo 



In 



(-) 



- 1 



AT 2 
12T 2 



(S15) 



The semi-phenomenological approximation quoted in Eq. (13) of the main text is obtained by dividing the second- 



order contribution by 2. As illustrated by the comparison with the full expression, Eq. (qT3j), in Fig. S2 the deviations 
of the approximate expression Eq. (13) from Eq. (£ 13) are insignificant at the present overall level of accuracy if AT 
is in the range . . . 150 K. 



Inserting Thbm from Eq. (£13) and the effective viscosity ?7hbm from Eq. (q6| into the generalized Stokes-Einstein 
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FIG. S2: Dependence of various expressions for the effecti ve Brownian temperature 7hbm on the temperature incre ment 
AT = T s ~ Tq: the calibrated scalar model given in Eq. (£131 (dot-dashed); the second order Taylor expansion Eq. (£151 
(dotted, blue); and the semi-phenomenological approximation given in Eq. (13), and displayed in Fig. 2, of the main text (solid, 
red). In contrast to the truncated Taylor series Eq. (£ |15[ ), the semi-phenomenological expression has no significant errors over 
the practically relevant temperature range. 



relation Dhbm = fcB^HBM/ (67t?7hbm-R)j Eq- (14) of the main text, we obtain the effective diffusion coefficient 

exp [-X(AT + T + T VF )/T VF ] 
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(S16) 



with abbreviations as above. As illustrated in Fig. |S3[ a satisfying approximation is indeed again obtained by use of 
the simple approximate expressions for Thbm and ijhbm from Eqs. (8) & (13) of the main text, respectively. This is 
how the curves in Fig. 2 of the main text were generated. 



4. APPARENT DIFFUSION COEFFICIENT FOR INHOMOGENEOUS HEATING POWER 



In practical applications diffraction usually limits the experimental realization of a spatially uniform heating rate 
within the observation volume. Therefore, the apparent diffusion coefficient deduced from the transit time statistics 
of the particles passing through the focus involves some implicit averaging over a spatially heterogeneous -Dhbm(f). 
In the following, we derive a theoretical expression for this average. 



First, the local diffusion coefficient -Dhbm(i") follows immediately from the heating power density /(r) via Eq. (£16) 
by noting that AT(r) a 7(r) if the absorption coefficient of the particle is temperature insensitive. Assuming a 
radially symmetric heating power distribution in the focus, the appropriate averaging procedure is similar to that 
for a particle released at the center of the focus, which can be traced back to a standard first-passage-time problem 
[S3! IS4| . The distribution P(r,t) of escape times for the particle is obtained by solving the Smoluchowski equation 
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FIG. S3: Temperature dependence of the effective diffusion coefficient Dhbm(AT) from the scalar model, as given in Eq. (S 161 
(dot-dashed), and in the approximation obtained by inserting the individual approximations of Thbm and ?7hbm from Eqs. (8) 
& (13) of the main text into the generalized Stokes-Einstein formula (solid, red). 



with an absorbing boundary condition at the boundary B of the focus volume, 

dP 

— =V-D(r)VP, P(r,0) =5(r-r ), P(r, f) = on B (S17) 

The spherically symmetric boundary value problem for the escape time r p (r) = r p (r) of a particle starting at position 
r is 

r p (2/r + d r )D + Dr^ = -1 with r p (r) = on B , (S18) 

which has the general solution 



^ ) = JAw^-mP)) ir '- (S19) 

uj being the radius of the focus volume and c a constant of integration. For the escape problem of a particle starting 
in the center of the focus, c = is required by t p (0) < oo. The apparent diffusion coefficient thus reads 

Dhbm = ^ 2 /6t- p (0) . (S20) 

For the related transit problem, which is of interest for our transit time analysis, the situation is slightly more 
complicated. The most likely transit paths are only touching or barely entering the focus, so that the transit time 
distribution Ptransit(^) diverges at t = 0. However, the characteristic transit time r t may be extracted from the 
experimentally obtained transit time distribution by fitting the asymptotic law Ptransit(^ S> 0) ~ i~ 3 / 2 exp[— i/r t ] 
[S51 IS6j . The stochastic errors inferred from the fits are displayed as error bars in Fig. 2 of the main text. The 
apparent diffusion coefficient follows from r t as 

^hbm = w 2 /^ 2 n , (S21) 

for spherical focus geometry. In practice, the focus is usually more elongated along the optical axis than transverse 
to it, so that it may be better approximated by a cylinder, in which case 

Dhbm = u 2 /ain , (S22) 
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where a% denotes the first zero of the Bessel function of the first kind of order zero. The different numerical factors 
in the last two expressions therefore provide a lower bound for the systematic numerical uncertainties involved in the 
determination of the absolute value of -Dhbm- 

A typical experiment yields a time series of photothermal bursts .ST. Their intensity is proportional to the 
local power densities of the lasers used for heating the particle and detecting the induced refractive index change, 
respectively. Here, an uncertainty arises since the focus geometry cannot be controlled or determined precisely in the 
diffusion experiment. A nominal lateral focus size p = 300 nm has therefore been estimated by fitting a Gaussian 
intensity profile cx exp[— r 2 /(2p 2 )] to the photothermal image of single immobilized gold nano-particles obtained with 
the same setup. As the axial extension of the focus is usually large compared to the lateral one (about 1 pm), the focal 



volume is approximated by a cylindrical shape, corresponding to Eq. ( 22 ) . Particles are identified as "in the focus" if 
the signal intensity surpasses a certain threshold set to a fixed percentage of the maximum signal of the whole time 
trace of bursts. The threshold therefore defines the actual focus size u relevant for the burst width analysis, which 
stays constant during the measurements of a given sample, due to scaling of the threshold with the maximum signal. 
In Fig. 2 of the main text, we use the value of u as a freely adjustable overall fit parameter to match the experimental 
data with the theoretical prediction for -Dhbm and find cj sa 250 nm for the R = 60 nm particles and oj s» 170 nm for 
the R — 40 nm gold nano-particles. Note that uj may generally differ between the measurements of different samples 
(viz. particle sizes) due to variations in the signal-to- noise ratio and the sample geometry. 

5. PARAMETERS FOR THE SOLVENT VISCOSITY 

Parameter values used for all graphs throughout the main text and the supplementary material: 

ambient temperature To = 298 K 

dynamic viscosity of water v{T) = Voo exp[A/ (T — Typ)] 

with r/oo = 0.0298376 mPa- s 

A = 496.889 K 
T VF = 152.0 K 
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